Collective Atomic Recoil Lasing Including Friction and Diffusion Effects 
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We extend the Collective Atomic Recoil Lasing (CARL) model including the effects of friction and 
diffusion forces acting on the atoms due to the presence of optical molasses fields. The results from 
this model are consistent with those from a recent experiment by Kruse et al. [Phys. Rev. Lett. 91, 
183601 (2003)]. In particular, we obtain a threshold condition above which collective backscattering 
occurs. Using a nonlinear analysis we show that the backscattered field and the bunching evolve 
to a steady-state, in contrast to the non- stationary behaviour of the standard CARL model. For a 
proper choice of the parameters, this steady-state can be superfluorescent. 



The mechanical effect of light on atoms has now been 
the subject of intense theoretical and experimental re- 
search efforts for several decades. However, the fact that 
the collective atomic centre-of-mass motion of atoms can 
strongly influence the evolution of optical fields has only 
received attention relatively recently 0, 0> El- Re- 
cent experimental studies involvinglarge numbers of cold 
atoms in high-quality cavities 0,00 represent an impor- 
tant advance in this field, allowing detailed experimen- 
tal studies of collective atom- light interaction dynamics. 
During these interactions both the mechanical effect of 
the cavity modes on the atomic motion and the driving 
of the cavity modes by the dynamic spatial distribution 
of atoms in the cavity must be described self-consistently 
and cannot be considered independently. 

A recent experiment by Kruse et al. represents the 
first unambiguous realization of the CARL model orig- 
inally proposed by Bonifacio and coworkers [jj, which 
describes collective backscattering of an optical pump 
field by a sample of cold atoms. Previous experiments 
on CARL have been performed in hot atomic vapours 
0, 0] where however the gain of the backward field can 
not be unambiguously attributed to atomic recoil. 

Here we extend the previous theoretical work on 
CARL, including the effects of friction and diffusion 
forces acting on the atoms due to the presence of op- 
tical molasses fields. We describe the system using a set 
of coupled Maxwell-Fokker-Planck equations. A linear 
stability analysis reveals that there is a threshold condi- 
tion for the pump power above which collective backscat- 
tering occurs. Preliminary experimental results confirm 
this prediction [lfj. We further show that the backscat- 
tered field and atomic density modulation amplitude or 
'bunching' evolve to a steady-state, in contrast to the 
non-stationary behaviour observed using the standard 
CARL model. Our model describes the main features 
of the experimental results of . 

In addition to the so-called 'good-cavity' regime in 
which the experiment of operates, we also examine 
the behaviour of the system in the 'bad-cavity' regime. 



We show that in this regime the atoms emit in a superflu- 
orescent way with scattered intensity cx N 2 , where 
N is the number of atoms. This is a unique example of 
a steady-state superfluorescence, i.e. superradiance from 
an incoherently prepared atomic system. 

Our system consists of an ensemble of atoms backscat- 
tering a far-detuned pump field into a counterpropagat- 
ing mode of a ring cavity, described by a set of equations 
derived by Bonifacio and coworkers 0, ■ The equa- 
tions are here generalized to include a friction force, — jp, 
and a stochastic force, F(t), due to the presence of optical 
molasses: 
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where i = u> r pt, 9 = 2kz, p = 2kv z /uj r p, and 
A = faoV/NhLup) 1 / 2 E are universally scaled time, 
position, momentum and scattered electric field vari- 
ables, respectively, normalized to the CARL bandwidth 
uj r p where uj r — 2fik 2 /m is the recoil frequency, p — 
(n /2A) 2 / 3 (u J d 2 N/2Vhe uj 2 r ) 1 / 3 oc (PqA/A 2 ) 1 / 3 , Q is 
the Rabi frequency of the pump field with frequency uj — 
ck, detuned from the atomic resonance by A = to — luq, 
d is the electric dipole moment of the atom, Po is the 
intracavity pump power and V is the cavity mode vol- 
ume. K — n c /uj r p represents scaled cavity losses and 
7 = "ff/Lu r p is the scaled damping coefficient to account 
for molasses friction. With the exception of the stochas- 
tic force, F(t), these equations are similar to the ones 
used by Bonifacio and Verkerk [l2( to describe CARL 
including the effect of collisions. 

Here, the stochastic force, F(t), causes diffusion of the 
atomic momenta i.e. heating. We assume white noise, 
i.e. (F(F)F(F)) = 2D p 5{t- F), where D p = ja 2 is the 
momentum diffusion coefficient and a is the momentum 
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spread in units of u r p, corresponding to the Doppler tem- 
perature T of the atoms in the molasses fields. 

In the limit of strong viscous damping, it is possible to 
adiabatically eliminate the atomic momentum for each 
atom by setting ^ = in Eq. J2J, so that 
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A simple steady-state analytical solution can be derived 
from Eqs.© and Q in absence of the stochastic force 
F, i.e. neglecting diffusion, and for the case of perfect 
bunching. In fact, assuming (e~ ie ) ~ e~ % ( 6 > in Eqs.|0 
and Q), one obtains p| 
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where we have normalized all the variables in order to re- 
duce the number of free parameters, defining a — A/y/j, 
p = y/jp = d9/dr, t — i/^/j and k = y/jK. In the fol- 
lowing we will use these scaled variables and parameters. 

Depending on the value of n, it is possible to define 
two different steady-state regimes for CARL: when n <C 
1 ('good-cavity' limit), we obtain \{p)\ w (2k) 1 / 3 ^> k 
and \a\ 2 sa (2k)" 2 / 3 . Since k cx k c ( 7/ A 2 /NP )^ 2 and 
| a, | 2 cx P s /N, where P s is the scattered light power, then 
P s cx N 4 / 3 as in the usual CARL 0,1111. Conversely, for 
/?> 1 ('bad-cavity' limit) we obtain \ (p)\ <C k and |a| 2 « 
1/k 2 , so that the scattered power P s is proportional to 
N 2 , i.e. superfluorescent. 

Diffusion can be described, in the adiabatic limit of 
Eq.Q, writing a Fokker-Planck equation for the distri- 
bution function P(9, t) [Hill. Together with Eq. ©, 
our model becomes: 
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where D = y/jD e = ° 2 / y/l oc T(A 2 /jfNPo) 1 / 2 and 
Dg = a 2 1 7 is the space diffusion coefficient. 

The distribution function is normalised such that 
J P{0, r) d6 = 1 and is periodic in 9 with pe- 
riod 2ir. Hence, Eq. JJjJl and Q can be writ- 
ten in terms of the spatial harmonics of P(9, r) i.e. 
P(9,r) = (1/27T) A, (rjeto*, whe re B n (r) = 

P(9,T)e~ m9 d9, so that 
dB n 



in (aB„_i + a*B n+1 ) - n 2 DB n (8) 
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We note that B_ n = i3* and Bq = 1. In particular 
& = |-Bx j = |(e~ j6l )| is the bunching factor, describing the 
ampitude of the density grating. Moreover, from Eq.Q it 
follows that the average momentum is (p) = — 2Re(a£?*). 



In the following, we will show that diffusion is responsi- 
ble for the existence of a definite threshold for the steady- 
state solution, for which in general the bunching factor b 
is less than one. 

The system of equations JSJ) and © have a steady- 
state solution with a*- -* = and Bn = for n ^ 0. 
If we linearise Eq. JHJl and © around this steady-state, 
we find that fluctuations of a and B n grow exponentially 
cx exp (At) if the dispersion relation 



(A + k) (A + D) = i 



(10) 



has roots with positive real parts. The real and imaginary 
parts of the unstable mode of Eq. (fTT)f) gives the gain G 
and the frequency shift Alu in the exponential regime, 
that, in units of the cavity bandwidth 
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where C = (k — D)/2. The threshold condition for insta- 
bility (G > 0) gives rise to the following relation between 
k and D: 



kD (D + k) < 1. 



(13) 



Figure shows the regions of parameter space (k,D) 
for which Eq. I|13|) predicts unstable growth of the probe 
field and the atomic density modulation. 




FIG. 1: Region of instability as a function of k and D (gray 
area) . 

From the experiment of Kruse et al. Q , the minimum 
value of the frequency shift measured was Auj/k c ~ 4.6. 
Assuming that it corresponds to the threshold value of 
the frequency shift, then, using Ea. H13l) in Eq. l|12|) we cal- 
culate that the values of k and D at threshold are k sa 0.1 
and D w 2.1. The experiment of Ref.Q can therefore be 
described using the 'good cavity' limit where k <C 1. In 
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this limit, the threshold condition Ea. H13|) becomes sim- 
ply D < k -1 / 3 , and the frequency shift at threshold is 
approximately Alo/k c « k~ 2 / 3 r* D 2 . Using the defini- 
tion of k, the threshold value k — 0.1 can be expressed 
as a threshold condition for thepump power. Assum- 
ing the same parameters as in p|, i.e k c = (27r)22kHz, 
7/ = 9k c and T w 150/i-ftT, we obtain a CARL parameter 
p = 14.6 at threshold, corresponding to an intracavity 
pump power of Pq m 3W, in good agreement with the 
experimental value. The experimentally accessible pa- 
rameters are the pump power and the CARL frequency 
at threshold. The linear analysis can provide us with 
important information on the scaling behaviour of the 
threshold input pump power and the frequency shift at 
threshold. Using the condition D < k -1 / 3 and the scal- 
ing of D and k given above, one obtains that the intra- 
cavity pump power scales as 



T 3 / 2 A 2 ^ /2 



(14) 



In the same 'good-cavity' limit, from Ea. (|12|) . the fre- 
quency shift at threshold can be shown to scale as 
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It is also worthwhile to investigate the scaling behav- 
ior of the threshold pump power and the frequency shift 
at threshold in the superfluorescent regime, k> 1. Al- 
though to date there have been no experimental studies 
which operate in this regime, we can use our model to 
predict the behaviour of such an experiment. It should be 
noted that from an inspection of the region of instability 
shown in figure^ hi order to operate in the superfluores- 
cent regime and remain above the threshold for instabil- 
ity, it is necessary to increase scaled cavity losses k and 
decrease the scaled diffusion/temperature parameter D. 
In the bad-cavity limit the threshold condition Ea. H13|) 
becomes simply D < k~ 3 , and the frequency shift at 
threshold is approximately Alo/k c « kT 2 w D 2 I Z . Con- 
sequently the threshold pump power in the superfluores- 
cent regime scales as 



T i / 2 A 2 «y 2 7 }/ 2 

N ' 



(16) 



Note that the dependence of threshold pump power on 
temperature is now cx T 1 / 2 in the bad cavity limit, as 
opposed to cx T 3 / 2 in the good cavity limit. The depen- 
dence on N and A is the same as that for the good cavity 
limit. The frequency shift at threshold in the bad cavity 
limit can be shown to scale as 
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so the dependence of the frequency shift at threshold on 
temperature, friction and cavity losses will be the same 
in both the good and bad cavity limits. 



The numerical integration of Eqs. J3J and © shows 
that the system evolves toward the stationary solution 
with at*' = ae luJT (where to = dtfi/dr is the scaled fre- 
quency shift) and = /3 n e muT ', where a and j3 n 
are complex constants. As a consequence, the den- 
sity distribution moves at a constant velocity, P{6, r) = 
p( s \9 + lot). We note that at steady-state the average 
momentum is (p) = — 2Kb 2 /(K 2 +co 2 ), where the bunching 
parameter is 6 = \/3i\. Comparing with Eq.JSJ, it follows 
that in the case of perfect bunching (6 = 1) to = — (p). 
Figure El shows the evolution of the backscattered scaled 
intensity, \a\ 2 , (a), and the bunching 6, (b), as a function 
of scaled time r for parameters close to those of the ex- 
periment by Kruse et al. i.e. D = 1.49 and K — 0.075. 
The instability was initiated using a seed field ao = 10~ 5 . 

It can be seen that the field intensity and the bunching 
increase exponentially before relaxing to a steady-state. 
This is the same qualitative behaviour as observed in 
the experiment 0- Fig[5] (c) shows also the time evo- 
lution of the frequency shift Aoj/k c = uj/k (continuous 
line) and of the scaled mean velocity 2k(v z ) / k c — (p) / n 
(dashed line). We observe that the atomic mean velocity 
at steady-state does not coincide with the velocity of the 
optical standing wave, i.e. — (p) =/= to, as occurs in the 
case of perfect bunching 6=1. Finally, Fig|2] (d) shows 
the stationary distribution P(0) vs. 8, showing the den- 
sity grating profile. The steady-state solution of Eqs. © 
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FIG. 2: Graph of |a| 2 , (a), b, (b), ui/k (continuous line) and 
— {p)/k (dotted line), (c), vs. r and of the stationary distri- 
bution P(6) vs. 9, (d). 

and may be obtained solving the following recurrence 
equation for (3 n : 



(u> - inD)/3 n = a/3„_i + a*(i n 



+1j 



where n ^ and a — /3i/(k + ilo). Calculating [3\ in 
terms of a continued fraction and iterating numerically 
the solution in order to find to, the stationary solution 
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can be obtained exactly. As an example, FigOt shows 
the steady-state bunching b (continuous line) and FiglJb 
shows —{p)/n (continuous line) and uj/k (dashed line) 
as a function of D and for n = 0.1. We observe that 
the bunching parameter b goes to zero at the threshold 
value Dth = 2.1, and the frequency shift at threshold 
is u = 4.6k, in agreement with the results of the linear 
theory. We note from Fig|3Jb) that — (p) w uj only well 
above threshold, when D <C D th- 
in this limit an approximate solution for the frequency 
shift uj and the bunching parameter b can be obtained 
assuming that the steady-state probability distribution 
P( s ) (9 + ujt) is a Gaussian. In this case it is possible to 
show that b and uj can be obtained from the solution of 
the following two coupled equations: 



2Kb 1 



UJ 



exp 



2V2u 



(18) 
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The approximate solution of Eas. (|18fl and (|19|) for b 
(dashed line) is shown by a dashed line in figEIa) for 
comparison with the exact solution. We note that the 
approximate solution is valid only when — (p) ~ uj, as 
can be seen from fig 03. 
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FIG. 3: (a): Steady-state solution for b (continuous line). 
The dashed line represents the approximated solution 1191 . 
(b): Steady-state solution for uj/k (dashed line) and —{p)/k 
(continuous line), as a function of D for k — 0.1. 

In the experiment of ref. 5] , the pump power is ramped 
up and down, while the frequency shift uj/k is monitored. 
The ramp can be designed in such a way that the pump 
power crosses the threshold value Pp, so that we expect 
that the signal of power P s of the scattered light beam 



has the dependencies shown in Fig. 0] This behavior 
seems to be in agreement with preliminary experimental 
results, to be discussed elsewhere [Io|. 

In conclusion, we have presented a model of the re- 
cent experiments by Kruse et al. Q showing collective 
backscattcring of an optical pump field by a sample of 
cold atoms. The model used is an adaptation of the 
CARL model which includes the effects of friction and 
diffusion forces acting on the atoms due to the presence 
of optical molasses fields. Using this model, we show 
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FIG. 4: Graph of frequency shift uj/k, (a), and scattered 
power | <x | 2 , (b), as a function of Pq/Pt, for k c — (2ir)22kHz, 
j f = 9k c and T = 150fiK. 



that the system can be described by a system of cou- 
pled Maxwell- Fokker-Planck equations. It was shown 
using a linear stability analysis that there is a thresh- 
old condition above which collective backscattering oc- 
curs. Using a nonlinear analysis it was shown that the 
backscattered field and atomic density modulation am- 
plitude or 'bunching' evolves to a steady-state, in con- 
trast to the non-stationary behaviour observed using the 
usual CARL model. We suggest to test our model ex- 
perimentally as well in the linear regime by searching for 
a threshold and characterizing its dependence from the 
coupling parameter and from temperature, as in the non- 
linear regime. The experiment currently, operates in the 
so-called 'good-cavity' regime We furthermore pro- 
pose to tune the experiment into the 'bad-cavity' regime 
in which the atoms scatter superfluorescently, with scat- 
tered intensity oc N 2 . This would represent the first cw 
superfluorescent system realized so far. 
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